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A GENERAL FORMULA FOR THE DECAY LIFETIMES 
OF THE STARFISH ELECTRONS* 


by 

E. G. StassinopoLilos 
Goddard Space Flight Center 
and 

P. Verzariu 

Applied Physics Laboratory, Johns Hopkins University 


INTRODUCTION 

Seven artificial radiation belts have been made by the explosion of high-altitude nuclear devices 
since 1958. Of these, the Starfish event of July 9, 1962, a low-latitude nuclear explosion, produced an 
intense region of trapped artificial electrons with very long lifetimes. 

The spatial and temporal characteristics of the artificial radiation zone have received considerable 
treatment in the literature, and the decay lifetimes for the injected electrons have been discussed by 
several authorsfVan Allen, 1 964; Cladis et ah, 1965; Bostrom and Williams, 1965;Brown, 1966; 
Mcllwain, 1 966; Paulikas et al., 1967; Beall et ah, 1967; Pfitzer, 1 968; Bostrom et al., 1970). The 
picture that emerges is that of a fast decay at very low L values (L ^ 1.2), which is consistent with the 
predictions of atmospheric scattering theory (Walt, 1964). Short lifetimes also occur at the outer edge 
of the inner radiation zone in the slot region. For the intermediate L values, the decay is slower, with 
longest lifetimes occurring at L ~ 1 .4. At low L values ( 1 . 1 5 ^ L ^ 1 .4), a ^-dependence is also ob- 
served (Beall et al., 1967). 

In the years following 1962, computer programs were developed at GSFC to estimate the tlux 
decay of artificially injected electrons and predict the radiation hazards to which satellites would be 
exposed (Stassinopoulos, 1965, 1967). The decay lifetimes used in the calculation of orbit-integrated, 
vehicle-encountered particle fluxes were based on fragmentary published data, and as such they were 
restricted and not very accurate (Stassinopoulos, 1965). In this paper a unified treatment of the 
decay of the artificial electrons in the inner zone is presented, with the decay lifetime r expressed as a 
continuous function of field strength B, magnetic shell parameter L, and energy E. The data used in 
the analysis is essentially taken from a single satellite, the 1963 38C,t and covers the time span from 


*This work was supported in part by the Department of the Navy under contract N00017-62-C-0604. 
t Designed and built by the Applied Physics Laboratory, Johns Hopkins University. 
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launch, September 28, 1963, to December 1968, when, it is felt, the energy channels of the satellite 
measured essentially the radiation of the natural background (Bostrom et al., 1970). A description of 
the satellite and its instrumentation can be found in the literature (Williams and Smith, 1965; Bostrom 
et al., 1967; Beall et al., 1967). 

DISCUSSION 

In order to facilitate the handling of the problem, the following simplifying assumptions are 
made: 

(1) The decay lifetime t remains constant in time. 

(2) Solar-cycle, diurnal, and secular-change effects can be neglected. 

(3) The decay of the Starfish electrons can be approximated by an exponential form. 

(4) The 5-dependence of r does not extend beyond L > 1 .40. 

Assumption (3) is justified by the 1963 38C data; assumption (4) is supported by comparison with 
OGO data, which measures a wider range of equatorial pitch angles than does 38C (Pfitzer, 1968; 
Pfitzer and Winckler, 1968). 

When the data is ordered in terms of the variable L, three distinct regions emerge (Figure 1), each 
identified by a characteristic functional dependence of r on the parameters B, L, and E. We denote 
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the three regions as I, II, and III. In region I, which extends from L= 1 . 1 5 to Z, = 1 .40 (for all electron 
energies), r is a function of all the variables; B, L, and E. The 5-dependence of the decay constant has 
been previously discussed by Beall et al. ( 1 967), and it appears to become indeterminate for L ^ 1 .4. 
Consequently, in our treatment we will take r in region II to depend only upon L and E. Region II 
extends from Z, > 1.4 to approximately L ~ 2.3, with the latter boundary being dependent on energy. 
Curves obtained for region II from the 38C data when t is plotted versus E for all three channels at 
given L values suggest that the high-energy electrons exhibit the most rapid decay and that at the time 
of the 1963 38C measurements almost no artificial electrons with energies greater than 2 MeV re- 
mained in regions where £^1.8. This is in general agreement with currently prevailing opinion* 
(Pfitzer, 1968). Thus, in the generation of the r versus E plots for region II (1 .4 < Z, ^ 2.3), no data 
from the 1963 38C channel 3 {E^ > 2.4 MeV) were used; instead, the r-dependence was extended lin- 
early to the higher energies, with only data obtained from channels 1 {E^ > 0.28 MeV) and 2 (£^ > 1 .2 
MeV) being used. This is consistent with the derivation of the decay lifetimes obtained by Beall et al. 
(1967) and reflects the uncertainty in the channel 3 data due to background proton contamination. 

At the higher L values (Z, > 1.8) in this region, the data were adjusted to fit reported measurements 
(Van Allen, 1 966; Mcllwain, 1966; Brown, 1966). 

Finally, in region III, which includes all Z-space beyond about 2.3, the decay lifetime was taken 
to be only a function of energy. In this region, r is typically of the order of one month or less and 
will be neglected in this treatment. Figure 1 outlines the domains of the regions discussed above for 
electrons with E^ > 0.28 MeV and indicates the dependence of t on the variables appropriate to each 
domain. 

In the formulation of this problem, a four-dimensional space is considered in which the decay 
lifetime r is a continuous function of the variables B, L, and E. The boundaries beyond which compu- 
tation of r is discontinued are: 

(1) The equatorial envelope, as defined by Mcllwain (1961). 

(2) The radiation-belt terminus, designated as the constant-electron-intensity contour of one 
particle per square centimeter per second, for trapped electrons with energies greater than 0.5 MeV, as 
given by Vette et al. (1966). 

(3) The surface t = 0. 

It is assumed that the terminus is independent of time and energy in 5Z,-space, since any slight 
changes in its position due to energy dependence or secular variations do not affect the derivations and 
conclusions presented in this treatment. The four-dimensional space is then subdivided into volume 
elements such that within each volume the parameter r may be linearized for each of its variables (i.e., 
all partial derivatives of r arc constant). 

Since the decay constant depends mainly on the magnetic shell parameter Z, appropriate incre- 
ments in Z have to be taken so as to preserve the linearity requirement for the derivatives of t within 
each volume element. Thus, region I is subdivided into 1 1 Z-intervals and region II into 19 Z-intervals. 


*J. I. Vette, private communication, 1970. 
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Figure 2— Decay lifetimes in region I (1 .15 < Z. < 1.4) showing the fi and Z. dependence of r for the 1963 38C satellite 

energy channels. 
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Figure 2 shows the dependence of t on the variables in region 1 for the three energy channels of satel- 
lite 1963 38C. In the next two figures, several cross sections of the t domain in region I are plotted. 
Figure 3 shows B-L cross sections for > 0.25, 0.5, 1 , and 3 MeV; Figure 4 shows several B-E cross 
sections at four discrete values of Z. In both figures, the dependence of r on B at these low values of 
L is clearly evident. Figure 5 shows the dependence of r on L and E in region II. The slowest decay in 
this treatment occurs at Z. ~ 1.6 for E^ ~ 0.3 MeV. The peak in the t contours gradually drifts toward 
T ~ 1.4 at the higher integral energies. Curves for r are not plotted for E^ < 0.2 MeV because there is 
evidence that the low-energy component (E^ ~ 0.1 MeV) of the integral electron flux does not mani- 
fest a perceptible decay (Pfitzer, 1968; Pfitzer and Winckler, 1968). 


DATA PROCESSING 

A generalized treatment for all three variables as they pertain to region I will be given here. The 
same treatment, appropriately modified to deal with only two variables, applies to region II. In general, 

t=t(B,L,E). (1) 

The following relations then hold within each volume element: 

T = k(B, L)E + MB, L), K = = ix{L)B + v{L), 


V = jL + T] , 
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L (earth radii) 

Figure 3~BL cross sections of the r domain in 




L (earth radii) 

I (1.15 <Z. < 1,4) for several integral energies. 


and 


A - + a{L) , a = eL + (p , 


9A , 


13 = 


dLdB ■ 


(3) 


Combining the expressions yields a general function for the mean lifetime within the volume element: 


t(B,L, E) — cxBLE + ^BL + "yLE + 5BE + eL + ^B + r\E + 0 . (. 

The coefficients of the variables B, L, and E form a matrix and are found for a particular volume ele- 
ment by successively linearizing all the derivatives of r. For those intermediate volume elements for 
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Figure 4— fiF cross sections of the t domain in region I ( 1 .15 < A < 1 .4) at several L values. 


which coefficients cannot be determined due to sparse or insufficient data points, the known coeffi- 
cient values from adjacent bracketing volume elements are linearly extrapolated, so as to match all 
slopes at the respective boundaries. Thus, a continuous functional relationship between the decay life- 
time T and its variables B, L, and E is generated for the inner zone. Table 1 gives the values of the 30 
sets of coefficients in the respective L ranges. In Table 2 the predicted and measured counting rates of 
channels 1 and 2 are compared at selected L values, for the period from January 1964 to December 1964. 
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Table 1 —Coefficients of the lifetime function. 


Volume 

Element 

Z,-Range 

a. 


7 

8 

e 


V 

0 


1 

1.15 -1.16 

6714.00 

24830.00 

-1236.00 

-8529.46 

-4439.40 

-31067.70 

1532.88 

5778.45 

2 

1.16 -1.17 

6322.00 

19470.00 

-1150.00 

-8074.70 

-2174.50 

-24849.90 

1433.12 

3151.15 

3 

1.17 -1.18 

8323.00 

15000.00 

-1579.00 

-10415.91 

-1500.00 

-19620.00 

1935.05 

2362.00 

4 

1.18 -1.20 

6584.50 

12080.00 

-1176.50 

-8364.48 

-963.30 

-16174.40 

1460.10 

1729.36 

5 

1.20 -1.25 

4203.60 

11854.00 

-776.20 

-5507.40 

-1070.74 

-15903.20 

979.74 

1858.29 

6 

1.25 -1.28 

3520.67 

10383.30 

-657.13 

-4653.73 

-878.33 

-14064.80 

830.91 

1617.78 

7 

1.28 -1.30 

2770.00 

9020.00 

-520.80 

-3692.88 

-863.25 

-12319.80 

656.40 

1598.48 

8 

1.30 -1.32 

1812.00 

8030.00 

-406.30 

-2447.48 

-820.80 

-11032.80 

507.55 

1543.29 

9 

1.32 -1.34 

1485.50 

7045.00 

-341.70 

-2016.50 

-506.20 

-9732.60 

422.28 

1128.02 

0 

1.34 -1.38 

578.00 

5634.15 

-139.25 

-800.45 

-484.33 

-7842.06 

151.00 

1098.71 

1 

1.38 -1.40 

140.50 

3346.70 

-18.00 

-196.70 

-216.95 

-4685.38 

-16.33 

729.73 

2 

1.40 -1.425 

0 

0 

-258.88 

0 

560.00 

0 

320.90 

-358.00 

3 

1.425-1.45 

0 

0 

-292.00 

0 

400.00 

0 

368.10 

-130.00 

4 

1.45 -1.475 

0 

0 

-276.00 

0 

320.00 

0 

344.90 

-14.00 

5 

1.475-1.51 

0 

0 

-297.14 

0 

200.00 

0 

376.09 

163.00 

6 

1.51 -1.55 

0 

0 

-337.75 

0 

175.00 

0 

437.40 

200.75 

7 

1.55 -1.60 

0 

0 

-422.22 

0 

160.00 

0 

568.33 

224.00 

8 

1.60 -1.625 

0 

0 

-425.20 

0 

40.00 

0 

573.10 

416.00 

9 

1.625-1.65 

0 

0 

-419.32 

0 

-80.00 

0 

563.55 

611.00 

10 

1.65 -1.70 

0 

0 

-358.34 

0 

-220.00 

0 

462.93 

842.00 

11 

1.70 -1.75 

0 

0 

-321.38 

0 

-400.00 

0 

400.10 

1148.00 

12 

1.75 -1.775 

0 

0 

-241.24 

0 

-520.00 

0 

259.85 

1358.00 

13 

1.775-1.82 

0 

0 

-147.78 

0 

-644.44 

0 

93.96 

1578.89 

14 

1.82 -1.85 

0 

0 

0 

0 

-766.67 

0 

-175.00 

1801.33 

15 

1.85 -1.87 

0 

0 

75.00 

0 

-850.00 

0 

-313.75 

1955.50 

16 

1.87 -1.95 

0 

0 

198.41 

0 

-962.50 

0 

-544.53 

2165.88 

17 

1.95 -2.00 

0 

0 

222.60 

0 

-980.00 

0 

-591.70 

2200.00 

18 

2.00 -2.05 

0 

0 

217.76 

0 

-900.00 

0 

-582.30 

2040.00 

19 

2.05 -2.10 

0 

0 

212.18 

0 

-740.00 

0 

-570.86 

1710.00 

30 

2.10 -2.30 

0 

0 

69.28 

0 

-600.00 

0 

-270.49 

1418.00 













CONCLUSIONS 


When the model presented in this treatment is used, the following should be considered: 

(1 ) Decay lifetimes should not be extrapolated backward too close to the epoch of the nuclear 
event because the initial decay (first few weeks) is considerably faster than at the later times considered 
in the present analysis. 

(2) The estimated maximum value of the error for the calculated t is about 2 percent at inter- 
mediate L values (1.2 ^ T 1.5). At higher L values (L ^ 1.6) and at low L values (L < 1.2), the 
error increases to about 30 percent on the average. 

(3) The functional dependence of r on energy has been derived through use of integral energy 
channels and may differ from decay constants obtained from differential energy measurements. 

(4) At low energies ^ 0.20 MeV), the validity of the energy dependence of t becomes mar- 
ginal. There is evidence that electrons of these energies do not manifest a decay (Pfitzer, 1968; Pfitzer 
and Winckler, 1968) but that they are mainly part of the natural component of the radiation belt. At 
the higher integral energies, it is felt that the treatment should not be extrapolated beyond > 5 MeV. 
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